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Annihilating dark matter (DM) has been discussed as a possible source of gamma-rays from the 
galactic center (GC) and contributing to the extragalactic gamma-ray background (EGB). Assuming 
universality of the density profile of DM halos, we show that it is quite unlikely that DM annihilation 
is a main constituent of EGB, without exceeding the observed gamma-ray flux from the GC. This 
argument becomes stronger when we include enhancement of the density profiles by supermassive 
black holes or baryon cooling. The presence of substructure may loosen the constraint, but only if 
a very large cross section as well as the rather flat profile are realized. 
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We have made great progress in our knowledge of what 
the universe is composed of. Surprisingly, we already 
know that the main constituents of the universe are not 
baryonic, but are unknown dark matter (DM) and dark 
energy. Recent analyses using observational data of the 
cosmic microwave background anisotropy, type la super¬ 
novae, and large scale structure precisely give the relic 
density of these components, Ddm = 0.22 and Da = 0.73 
[l]. In particular for DM, we have some candidates mo¬ 
tivated by particle physics. The most viable candidate 
is the lightest supersymmetric particle, which is a neu- 
tralino in most models. This supersymmetric neutralino 
can annihilate into final states including photons via var¬ 
ious channels, and these photons might be detectable or 
might have already been detected from several astrophys- 
ical objects (see Refs. [3] for reviews). 

In the direction of the galactic center (GC), there are 
clear gamma-ray signals in the GeV and TeV energy 
regions, which have been detected respectively by the 
EGRET detector Q, and the atmospheric Cerenkov tele¬ 
scopes (ACTs) such as Whipple @,0|, CANGAROO-II 
{fill, and H.E.S.S. .7]- These gamma-rays from the GC are 
potentially due to DM annihilation, and have been exten¬ 
sively studied Js) Q. These results show that if the den¬ 
sity profile of the galactic central region is cuspy enough, 
as suggested by IV-body simulations such as by Navarro, 
Frenk, and White 0 (hereafter NFW) and Moore et 
al. ljj (hereafter M99), then the gamnra-ray fluxes can 
be explained by the neutralino annihilation with a cross 
section that gives the proper relic density Ddm- 

On the other hand, analyses of the diffuse EGRET 
emission shows the signature of an extragalactic gamma- 
ray background (EGB) in the GeV range [12, |jjj| ■ An¬ 
nihilating DM may also significantly contribute to this 
EGB flux. Using the hierarchical clustering formalism 
that is now widely accepted, several authors gave the 
flux predictions, investi gati ng the effect of DM cluster¬ 
ing or substructure 0 , 0 , [la, Il3 , and suggested that 
the EGB data can be explained well by including the DM 
component. In particular, a bump around 3 GeV, discov¬ 
ered by the recent reanalysis yjj, may be the signature 
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In this Letter, we investigate the gamma-ray signa¬ 
ture from the GC and EGB together, assuming that the 
halo profile is universal as suggested by recent V-body 
simulations 0 - With this self-contained treatment, we 
point out that the annihilating DM cannot be a main con¬ 
stituent of the observed EGB without exceeding observa¬ 
tional bounds imposed by gamma-ray measurements of 
the GC. Since both the GC and EGB fluxes should be 
predicted using the same cross section and mass of the 
DM particle, they are connected if we specify these in¬ 
gredients. We also show that the main conclusion of this 
Letter has a quite robust characteristic such that it does 
not depend on uncertainties concerning both the particle 
physics models and other astrophysical inputs. The lat¬ 
ter includes the central spike of halos due to presence of 
supermassive black holes (SMBHs) 0 or to the baryon 
cooling [0, and the significant enhancement of the EGB 
flux due to the inclusion of halo substructure. 

Gamma-rays from the galactic center.- The number 
flux of high-energy gamma-rays due to DM annihilation 
can be calculated with the following formulation: 

$° C (£ 7 )AD = 9.4 x 1CT 11 cm -2 s" 1 to“ 2 2 

x <au)_ 26 ^J(AD)AD, (1) 
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where m x = 100m Xi 2 GeV is the mass of the DM par¬ 
ticle, (av) = 10 _26 (cru)_26 cm 3 s _1 is average value of 
the annihilation cross section times the relative velocity 
(assumed to be independent of v), and dN 1 /dE 1 rep¬ 
resents the gamma-ray spectrum per annihilation, for 
which we use a simple parameterization as dN~JdE 7 ~ 
(0.73/m x )e -7 ' 76 ' E '’'/ m */p 7 /m x ) 1,5 + 0.00014] 0. (Al¬ 
though this parameterization may be less precise, it is 
sufficient for our purpose.) J(AD) represents the aver¬ 
age value of the following quantity: 
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over the detector angular resolution AD [AD ~ 2 x 
10~~ 3 (4 x 10 -5 ) for EGRET (ACTs)], and normalized to 
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TABLE I: Angular acceptance J(ASl) of the GC gamma-rays, 
and local values of enhancement factor /(0) for the EGB flux. 


Model J(2 x 10- 3 ) J(4 x 10~ 5 ) f(z = 0) 

NFW 7 x 10 2 5 x 10 3 2 x 10 4 

M99 6 x 10 4 4 x 10 6 2 x 10 5 


the local value. The integration is performed along the 
line of sight (l.o.s.) labeled by angle deviation d from 
the GC (t/j = 0 for the direction to the GC), and r is the 
direction to the integrated point from the GC. 

For the density profile of the DM halos p(r), we use 
the NFW and M99 models, which are characterized by 
the central slopes of 7 = 1 and 1.5, respectively, where 7 
is defined by p(r) oc r -7 for small radii. While the most 
recent iV-body simulations suggest there is no asymp¬ 
totic slope and do not give such a cuspy profile as M99 
(NFW may still be marginally consistent) , we use 
these two profiles as our reference models, in order to in¬ 
vestigate how our conclusion changes with the selected 
profile. Furthermore, considering some other physical 
processes, it would still be possible to obtain significant 
enhancement of the central density. The proposed mech¬ 
anism giving such a “spike” that may be steeper than 
the M99 profile is the accretion of DM particles onto a 
central SMBH [l9] or the effect of baryon cooling [i(J ■ 

In Table |U we summarize the values of J(AD) evalu¬ 
ated with the NFW and M99 density profiles. Because of 
its steeper profile in the central region, the M99 profile 
gives much larger values of J(AS2) than NFW. This dif¬ 
ference becomes more prominent when we use detectors 
with better resolution, since a more concentrated region 
can be probed. In Table |T| we used the cutoff scale 10 -8 
kpc, and the profile was assumed to be flat within that 
radius. This is because without this implementation, the 
l.o.s. integration would diverge mathematically at very 
small radius. The cutoff scale may be physically deter¬ 
mined by the annihilation itself, scattering of DM particle 
off stars, contraction of baryons, or the presence of a cen¬ 
tral SMBH. Because even the most recent simulations do 
not reach the very inner region of the halo (but as large as 
~ 0.1 1 kpc), the choice of the cutoff scale is a nontrivial 
problem. However, the GC flux is rather weakly depen¬ 
dent on this parameter in the reasonable range [tj, and 
the uncertainty does not strongly affect our conclusion. 

Figure Cta) shows the gamma-ray flux from the GC 
due to annihilating DM, with masses 50 GeV or 2 TeV. 
In deriving these expressions, we assumed (av }-26 = 3 
that is considered to be appropriate for leaving the ob¬ 
served relic density of DM Q; larger values than this 
would imply a lower relic density, requiring an additional 
DM component. Data points in 0.03-10 GeV and 0.2-2 
TeV are taken from the EGRET and CANGAROO-II 
{g] papers. A solid line above 2 TeV represents the power- 
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FIG. 1: (a) Gamma-ray flux from the GC from annihilating 
DM, with mass 50 GeV or 2 TeV, evaluated with the NFW 
and M99 profiles. Data from EGRET 13] and CANGAROO- 
II (| are also plotted; the H.E.S.S. result Q is shown as a 
solid line, (b) EGB intensity from DM annihilation. EGRET 
data points 1131 are also plotted. 


law fit to the H.E.S.S. data [?f (the recent Whipple result 
is consistent with the H.E.S.S. data 0 )- Correspond¬ 
ingly, theoretical curves are evaluated using Eq. 0 with 
A Q = 2 x 10 ~ 3 for 50 GeV, and with AO = 4 x 10 -5 for 
2 TeV DM particles. We can clearly see that, in the case 
of m x = 50 GeV, the flux evaluated with the M99 pro¬ 
file can easily be quite consistent with the EGRET data 
points over the wide range of energy. With the NFW 
profile, on the other hand, we predict considerably less 
flux. TeV gamma-rays may also be dominated by DM 
component, although both the flux and spectral shape 
are still controversial. 

Extragalactic gamma-ray background .—The EGB flux 
estimation involves somewhat more information, e.g., 
that on the cosmological clustering of DM halos. The 
intensity of EGB is calculated by 
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where = (1 + z)E 1 , h{z ) = [(1 + z ) 3 + Da ] 1 / 2 
and perit is the critical density. We also include the ef¬ 
fect of gamma-ray absorption by e _r , which is caused 
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by pair annihilation with the diffuse extragalactic back¬ 
ground light in the infrared or optical wavebands [2ill . 
This effect changes the EGB flux at TeV regions, but 
its extent is too small to affect the results. Hierarchi¬ 
cal clustering of the DM halos is included in an inten¬ 
sity multiplier f(z) 130. For evaluating it we used 
the halo mass function based on the ellipsoidal collapse 
model m with a lower mass cutoff of M m ; n = 10 6 Mq, 
which may be determined by the validity of hierarchi¬ 
cal clustering formalism, self-limitation due to annihila¬ 
tion itself, or nuclear and star formation activities [lfl . 
Varying it over a reasonable range (10 4 -10 8 Mm) changes 
the EGB flux only by a factor of 2 or less [0 12]. To 
evaluate a concentration parameter that represents how 
the bulk of mass in each halo concentrates in the cen¬ 
tral region, we used a publicly available numerical code 
by Eke, Navarro, and Steinmetz (23]. The resulting val¬ 
ues of /(0) for each profile are summarized in the fourth 
column of Table UJ and we note that our result is con¬ 
sistent with that of Refs. 00. We should note that 
these values are significantly smaller than those adopted 
by previous studies such as Refs. 0 . This large discrep¬ 
ancy potentially comes from uncertainty concerning the 
concentration parameter, and the presence or absence of 
substructures. We also discuss these possibilities later. 

We show in Fig. mb) the EGB intensity, with the same 
physical inputs as in the GC flux calculation. This shows 
that without any processes that give much larger f(z), 
the expected DM contribution to the EGB flux is consid¬ 
erably smaller than the observed value. We also note that 
the dependence on the adopted profile is less prominent 
in the case of EGB, compared with strong dependence of 
gamma-ray flux from the GC. This is because the EGB 
flux is less sensitive to the very central region of the halo. 
The weaker dependence of /(0) on the profile shown in 
Table d also reflects the same characteristic. With our 
canonical model, it is much more difficult to explain the 
observed EGB intensity mainly by annihilating DM com¬ 
ponent than gamma-ray flux from the GC; it requires an 
additional boost by more than two orders of magnitude. 

Constraints on annihilating dark matter component .— 
The contribution of DM annihilation to the EGB flux is 
quite strongly constrained by the GC gamma-ray obser¬ 
vations, and this result is rather robust, independent of 
uncertainties in the particle physics models. As a first 
step, we introduce a boost factor b for both the GC 
gamma-rays and EGB, as a correction to the canonical 
predictions of each flux. All the corrections due to the 
other astrophysical and particle physical possibilities are 
included in b. In order for the DM annihilation to be 
a main constituent of the observed fluxes, the required 
values of b should be very close to the following quan¬ 
tity: 6 max = miiij [$°^ > f/$^ h (£ : 7 ,i)], where i represents 
bin-number of each observation, and <f>° b f and are 
the observed intensity in ith bin and theoretical predic¬ 
tion given by Eqs. © and ®, respectively. By taking 
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FIG. 2: Boost factors, &Gcf an< ! ^egbi required to make DM 
annihilation a main component of the data, for the NFW and 
M99 profiles. Assumed masses are 10, 20, 50, 100, 200, 500, 
1000, 2000, and 5000 GeV (from bottom to top). The solid 
and dotted lines represent 6egb = &gc, 10&gc, 100&gc and 
&GC = 1, 6egb = 1, plotted for comparison. 

minimum over all the bins, we renormalize the flux with 
keeping its shape, so as not to exceed the data points, 
which should be regarded as rigorous upper limits. 

From Fig. ^ we can confirm that the value of b re¬ 
quired for the GC gamma-ray data is much smaller than 
that for EGB, i.e., 6 qq x <C &egb' These rec ffured values 
(6 qq x , ^eqe) are plotted in Fig. |3 for both the NFW and 
M99 profiles and for several assumed DM particle masses. 
We used the CANGAROO-II data in the TeV region, 
since this would be more conservative for our purpose. 
If we add corrections that are related only to particle 
physics models (especially by renormalizing the cross sec¬ 
tion), we obtain the relation between two boost factors as 
6egb = 6gc, since the correction is common to the both 
cases. In this case, because of the relation 6 egb ^ ^Gcf 
for all the models as shown in Fig. [3 we cannot explain 
the EGB data mainly by annihilating DM. Otherwise, it 
would overproduce gamma-rays from the GC compared 
to the data. In addition, Fig. |3also shows that this ten¬ 
dency is more prominent for the M99 profile than NFW. 
Although we restrict our argument within these two spe¬ 
cific profiles, the conclusion derived here is general and 
applicable to other profile choices, as discussed below for 
more specific examples. 

We note that the strong gamma-ray signal may not 
be coming from the GC; it has been suggested that the 
EGRET GeV source position is offset from the dynamical 
center of the galaxy, Sgr A*, at roughly 95% C.L. [0. 
If this is true, it also strengthens our argument, because 
it suggests that the most of the gamma-rays come from 
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other astrophysical sources, and DM component should 
be significantly smaller than the EGRET data. 

Other astrophysical possibilities.- A SMBH, due to its 
deep potential well, could accrete a significant amount 
of DM particles, and this would make the density spike 
in the central region of halos 19]. It has also been 
pointed out that the infall of baryons due to radiative 
cooling could lead the DM compression in the GC |2flj . 
Both these effects, potentially and significantly, enhance 
gamma-ray signals from DM halos. It should be noted 
that an enhancement of the central density profile by any 
possible processes strengthens our main conclusion. This 
is because the GC gamma-ray flux is much more sensitive 
to the slope of the central region, and the resulting rela¬ 
tion of the boost factors, bcc > ^egb, prevents the DM 
component from becoming dominant in EGB, without 
violating the gamma-ray observations of the GC. 

On the other hand, if the density profile in the central 
region of halos is less steep than the NFW (due to, e.g., 
rather large inner cutoff radius, as already mentioned), 
the required relation between the boost factors would 
become close to &egb ~ ^Gcf- I n this case, however, 
we should note that the DM component becomes abso¬ 
lutely difficult to be dominant both in the GC and EGB, 
while it is relatively easier to explain the EGB flux with¬ 
out overproducing the GC gamma-rays; it requires, e.g., 
much larger cross section, which is unlikely. For exam¬ 
ple, the calculation using the profile with 7 = 0.5 and 
rn x = 100 GeV gives b ~ 56gg x = 3 x 10 3 . 

Recent iV-body simulations suggest the presence of 
DM substructure, although it is not observationally con¬ 
firmed. According to the literature, this might boost the 
GC gam ma-ray and EGB flux by at most a factor of a 
few fr]| and about an order of magnitude 0 Q , respec¬ 
tively. Therefore, we obtain the relation, 6egb 1$ 106gc 
that is still below the required points shown in Fig. Q It 
suggests that even the inclusion of substructure cannot 
provide a way that violates the main thrust in this Letter. 
In the previous studies of the EGB flux E3 , the intensity 
multiplier as large as /(0) ~ 10' (for the M99 profile) was 
used, which is about a factor of 50 larger than our value 
(see Table QJ. This discrepancy may come from the dif¬ 
ferent choice of the concentration parameter, in addition 
to inclusion of substructure. The former is extensively 
discussed in Ref. 0 , and found to give an uncertainty 
of a factor ~ 5. Even if we use this extreme boost factor 
for the EGB (6egb ~ 505gc), for the M99 profile it still 
requires some additional effects that enhances the EGB 
flux with changing the GC gamma-rays by a significantly 
smaller amount. For the NFW or less steep profile, while 
it might provide a solution to the relative smallness of the 
predicted EGB flux, we still require considerable amount 
of corrections, which is physically unlikely. 
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